Chandra observations of the accretion-driven millisecond X-ray 
pulsars XTE J0929-314 and XTE J1751-305 in quiescence 

Rudy Wijnands^, Jeroen Homan^, Craig O. Heinke^'"^'^ Jon M. Miller^'^, Walter H. G. 

Lewin^ 

O : ABSTRACT 

o 

Oh' We observed the accretion- driven millisecond X-ray pulsars XTE J0929-314 

C/^ ■ and XTE J1751-305 in their quiescent states using Chandra. From XTE J0929- 

Q\ • 314 we detected 22 source photons (in the energy range 0.3-8 keV) in ~24.4 ksec, 

^^ ■ resulting in a background-corrected time-averaged count rate of 9±2 x 10~^ counts 

CN ■ s~^. The small number of photons detected did not allow for a detailed spectral 

j^^ ■ analysis of the quiescent spectrum, but we can demonstrate that the spectrum is 

^j^ • harder than simple thermal emission which is what is usually presumed to arise 

^ ■ from a cooling neutron star that has been heated during the outbursts. Assuming 

2 ' a power-law model for the time-averaged (averaged over the whole observation) 

O ■ X-ray spectrum, we obtain a power-law index of I.S^qj and an unabsorbed X-ray 

^ ■ flux of G'tt X 10~^^ ergs s~^ cm~^ (for the energy range 0.5-10 keV), resulting 
in a 0.5-10 keV X-ray luminosity of Ttl x 10^^ (rf/10 kpc)^ ergs s"^, with d the 



Dh' 



^ ■ distance toward the source in kpc. No thermal component could be detected; 

c/3 . such a component contributed at most 30% to the 0.5-10 keV flux. Variability 

in the count rate of XTE J0929-314 was observed at the 95% confidence level. 



% 



We did not conclusively detect XTE J1751-305 in our ~43 ksec observation, 
with 0.5-10 keV flux upper limits between 0.2 and 2.7 xlO^^^ ergs s~^ cm~^ 
depending on assumed spectral shape, resulting in 0.5-10 keV luminosity upper 
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limits of 0.2 - 2 xlO^^ ((i/8 kpc)^ ergs s~^. We compare our results with those 
obtained for other neutron-star X-ray transients in their quiescent state, and 
in particular with the quiescent properties of SAX J1808. 4-3658. Using simple 
accretion disk physics in combination with our measured quiescent luminosity 
of XTE J0929-314 and the luminosity upper limits of XTE J1751-305, and the 
known spin frequency of the neutron stars, we could constrain the magnetic field 
of the neutron stars in XTE J0929-314 and XTE J1751-305 to be less than 
3 X 10^YoT~c ^^^ 3 — 7 X lO^g-^-^ Gauss (depending on assumed spectral shape 
of the quiescent spectrum), respectively. 

Subject headings: accretion, accretion disks — stars: neutron stars: individual 
(XTE J0929-314; XTE J1751-305)— X-rays: stars 



1. Introduction 

Neutron stars in low-mass X-ray binaries accrete matter from solar mass companion 
stars. The sub-group of neutron-star transients spend most of their time in quiescence during 
which hardly any or no accretion occurs onto their neutron stars. However, these transients 
sporadically become very luminous in X-rays (>10^^~^^ ergs s~^) owing to a huge increase 
in the accretion rate. During those outbursts, these sources can be readily studied with the 
available X-ray instruments, but so far only about a dozen sources have been studied in their 
much dimmer quiescent states. In quiescence, they typically exhibit 0.5-10 keV luminosities 
of 10^^"^^ ergs s^^ and their spectra are usually dominated by a soft component which can 
be described by a thermal model (either a black-body or a neutron-star atmosphere model). 
This emission is generally ascribed to the cooling of the neutron star which has been heated 
during the outbursts (e.g., via deep crustal heating; Brown, Bildsten, & Rutledge 1998). 
For several quiescent systems an additional power-law shaped component is present in their 
X-ray spectra and it dominates above a few keV (e.g., Asai et al. 1998; Rutledge et al. 2001). 
The origin of this power-law component is not understood but it has been proposed to be 
due to residual accretion on the neutron-star magnetic field, an active pulsar mechanism, 
or shock emission due to the interaction of the pulsar wind and matter which is still being 
transferred from the companion star (e.g., Stella et al. 1994; Campana et al. 1998; Campana 
& Stella 2000). 

The fractional contribution of this power-law component to the flux in the 0.5-10 keV 
energy range varies significantly between systems. In some systems this component cannot 
be detected (with less than 10% of the 0.5-10 keV flux possibly due to such a component), 
but in other systems it contributes up to half the emission in the 0.5-10 keV energy range 
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(e.g., Asai et al. 1998; Rutledge et al. 2001). Currently, only two quiescent systems (SAX 
J 1808. 4-3658 and EXO 1745-248) have been found to exhibit quiescent spectra which are 
fully dominated by the power-law component (with over 90% of the 0.5-10 keV flux due 
to the power-law component; Campana et al. 2002; Wijnands et al. 2004). For several 
other systems (i.e., SAX J1810.8-2609 and XTE J2123-058; Jonker, Wijnands, & van der 
Klis 2004a; Tomsick et al. 2004) it has also been found that the quiescent spectra could be 
adequately fltted with only a power-law model. However, the statistics of the data for those 
sources were rather limited, still allowing for a thermal component which contributed more 
than 50%-70% to the 0.5-10 keV flux. 

During outburst SAX J1808. 4-3658 is an accretion-driven millisecond X-ray pulsar and 
its anomalous X-ray properties in quiescence (its hard quiescent spectrum and low luminosity 
of ~5 X 10'^^ ergs s^^; Campana et al. 2002) might be related to the expected higher magnetic 
fleld strength of the neutron star in this system, compared to that of the neutron stars in 
the non-pulsating systems. However, the recent study of the neutron-star X-ray transient 
EXO 1745-248 located in the globular cluster Terzan 5 during its quiescent state has cast 
doubt on this hypothesis. Wijnands et al. (2004) found that the Terzan 5 system, which 
is not observed to pulsate during outburst, also had a quiescent spectrum which was fully 
dominated by the power-law component (more than 90% of the 0.5-10 keV flux was due 
to the power-law component) just like SAX J1808. 4-3658. However, EXO 1745-248 had a 
quiescent luminosity (~2 x 10^^ ergs s'^) which was a factor of ~40 higher than what has 
been observed for SAX J1808.4-3658. 

The physical process(es) behind the power-law component and the differences among 
sources in terms of the properties of this spectral component are unknown. However, some 
clues about the nature of this spectral component were recently glimpsed in the work by 
Jonker et al. (2004a, b). They compared all quiescent neutron star X-ray transients for 
which spectral information was available and found that the fractional contribution of the 
power-law component to the 0.5-10 keV fluxes is lowest when these sources have quiescent 
luminosities of ~ 1 — 2 x 10^^ ergs s~^. They found that both at higher and lower quiescent 
luminosities the fractional contribution of the power-law component to the 0.5-10 keV flux 
increases. So far, only EXO 1745-248 did not follow this correlation. It remains to be seen 
whether this transient is an unusual system or if the correlation found by Jonker et al. (2004a, 
b) is spurious and more sources similar to EXO 1745-248 will be found. 

In recent years, four additional accreting millisecond X-ray pulsars (after SAX J1808.4- 
3658) have been discovered (Markwardt et al. 2002; Galloway et al. 2002; Markwardt, Smith, 
& Swank 2003; Markwardt & Swank 2003). Those systems are prime targets to observe in 
quiescence to test the hypothesis that the unusual quiescent properties of SAX J1808. 4-3658 
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are due to a relatively strong neutron-star magnetic field strength compared to that of the 
non-pulsating systems: one would expect that the magnetic fields of the neutron stars in 
those additional systems are of similar strength as that of the neutron star in SAX J1808.4- 
3658. To this end we had observations of XTE J1751-305 and XTE J0929-314 scheduled 
during cycle 5 of Chandra. Here we report on the results of those observations. 

The second accretion driven millisecond X-ray pulsar, XTE Jl 75 1-305, was discovered 
on April 3, 2002 (Markwardt et al. 2002). It was found to harbor a neutron star with a 
spin frequency of 435 Hz which is located in a binary system with an orbital period of 42 
minutes (Markwardt et al. 2002). The third accreting millisecond X-ray pulsar, XTE J0929- 
314, was discovered by Remillard (2002) on April 13-18, 2002, using the all-sky monitor 
(ASM) aboard the Rossi X-ray Timing Explorer [RXTE). Only later (May 2, 2002), when 
the proportional counter array aboard i^XTE' observed the source, was it discovered to be an 
accreting millisecond X-ray pulsar with a pulse frequency of 185 Hz (Remillard, Swank, & 
Strohmayer 2002). The orbital period of the system was found to be 44 minutes (Galloway 
et al. 2002). 



2. Observations, analysis, and results 

We observed XTE J0929-314 with Chandra on March 18, 2004, between 02:30 and 
10:02 UT (resulting in an exposure time of ~24.4 ksec) and XTE J1751-305 on June 26, 
2004 between 07:59 and 21:12 UT (resulting in an exposure time of ~43.0 ksec) We used 
the ACIS-S3 CCD and we selected a 1/4 sub-array to limit the pile- up in case the source 
fluxes exceeded ~10^^^ ergs s~^ cm^^ (as we will show below, both sources had quiescent 
fluxes significantly lower than this flux level so no pile- up occurred during our observations). 
We checked for background flares during our observations, but none were found allowing us 
to use all available data. Our analyzes were performed using the CIAO software package 
(version 3.1) and the standard Chandra analysis threads^. 



^See http://cxc.harvard.edu/ciao/ for CIAO and the analysis threads 
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2.1. XTE J0929-314 

2.1.1. Image analysis 

We used the tool 'wavdetect' to search for point sources in our data of XTE J0929-314 
and to obtain the coordinates of each source which was detected. For XTE J0929-314, we 
used the tool 'dmextract' to extract the number of observed photons and the count rate 
(for the energy range 0.3-8 keV). As source extraction region we used a circle centered on 
the source position (as obtained with 'wavdetect') with a radius of 1.5". For background 
extraction region we used an annulus centered on the source position with an inner radius 
of 5" and an outer radius of 20" (a larger outer radius could not be used since the field 
source CXOU J092919. 1-312244 was within 25" of XTE J0929-314 as can be seen in Fig. 1 
left panel). We detected four sources in our image, including XTE J0929-314; we will not 
discuss the other three sources further in this article. Instead we concentrate on the observed 
X-ray properties of XTE J0929-314 in its quiescent state. The coordinates we obtained for 
XTE J0929-314 are right ascension 09^^ 29™ 20M80 and declination -31° 23' 03".5 (epoch 
J2000.0), which are consistent with the coordinates of the source obtained during outbursts 
in the X-ray (Juett, Galloway, & Chakrabarty 2003), optical (Greenhill, Giles, & Hill 2002), 
and radio bands (Rupen, Dhawan, & Mioduszewski 2002; see Fig. 1 right panel). The total 
number of source photons detected in the energy range 0.3-8 keV from XTE J0929-314 is 
22, resulting in a net time-averaged count rate of 9 ± 2 x 10^^ counts s^^ (after background 
subtraction; less than 0.4 background photons are expected in our source extraction region). 



2.1.2. Count rate analysis 

In Figure 2 we plot the energies of the X-ray photons (0.3-8 keV energy range) detected 
at the source position against the arrival time of the photons (measured since the start of the 
observation). From this figure it can be seen that 4 photons are detected with energies above 
3 keV (note that no photons were detected with energies >5 keV) and the remainder of the 
photons have energies below ~2.5 keV. In the last ~10 ksec of the observation only two source 
photons are detected compared with twenty photons during the first ~15 ksec; this might 
suggest that the source exhibited variability during our observation. To investigate if XTE 
J0929-314 was indeed variable during our observation, we applied Kolmogorov-Smirnov and 
Cramer- Von Mises tests on the event list (as shown in Figure 2) to attempt to disprove the 
hypothesis that the source count rate is constant. Both tests showed that XTE J0929-314 is 
variable for energies between 0.3 and 8 keV (although effectively only up to 5 keV since no 
photons are detected above that energy) at the 95% confidence level. From Figure 2 it can 
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be seen that the four photons which have energies >3 keV all arrived within the first 5000 
seconds of the observation, which may be the main cause of the variability detected in the 
source. Therefore, we also applied the above two tests to the event list for photon energies 
between 0.3 and 2.5 keV. The Kolmogorov-Smirnov test revealed that the 0.3-2.5 keV event 
list might also be variable but only at a 90% confidence level. The Cramer- Von Mises test 
could not find evidence for variability in this energy range. Therefore, we conclude that the 
variability observed in the source is likely stronger at higher photon energies. 



2.1.3. Time- averaged X-ray spectrum 

Despite the low number of photons detected, we extracted the source spectra (averaged 
over the whole duration of the observation) using the CIAO tool 'psextract' for the energy 
range 0.3-8 keV, which also created the response matrix and the ancillary response files (the 
latter was also automatically corrected for the time-variable low-energy quantum efficiency 
degradation of the CCD). We used a circle with a radius of 1.5" as source extraction region. 
We fitted the spectrum using Xspec version 11.3.0 (Arnaud 1996). The small number of 
photons observed does not allow for x^ statistics to be used during the fits. Therefore, we 
fitted the data using Cash statistics (Cash 1979). Since Cash statistics cannot be used on 
background subtracted spectra, we did not subtract the background from our data. The 
errors thus introduced are likely to be small since less than 0.4 background photons are 
expected in the source extraction region (see above). The quahty of the fits was investigated 
by generating 10,000 Monte Carlo simulations of the best-fit spectrum: if the fit is good, 
roughly half the simulations should have values of the Cash statistics which are lower than 
those of the data. 

In all our spectral fits we left the interstellar column density A^h as a free parameter. 
To calculate the errors on the obtained fiuxes, we fixed each free fit parameter one at a time, 
either to its minimum or maximum allowed value as obtained from the fits. After that we 
refitted the data and recalculated the fiuxes. This process was then repeated for each free 
parameter and the final fiux range determined the fiux errors. The spectra of neutron-star X- 
ray transients in their quiescent states are usually dominated by a soft thermal component 
which can adequately be fitted with a neutron-star hydrogen atmosphere model for non- 
magnetized stars (the NSA model; we used that of Zavlin, Pavlov, & Shibanov 1996; we 
also assumed a mass of 1.4 Mq and a radius of 10 km for the neutron star in XTE J0929- 
314). In such a model the normalization is given by ^, with d the distance to the source in 
parsecs. The distance toward XTE J0929-314 is not known but it is constrained to be >5 kpc 
(Galloway et al. 2002). When leaving the normalization free, 88% of the simulated spectra 
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have better Cash-statistics than the data demonstrating that this model does not provide an 
accurate description of the data (Fig. 3 bottom). From this model we obtained an effective 
temperature kT^o (for an observer at infinity) of 0.3 ±0.1 keV and the column density is 
<3 X 10^*^ cm^^. However, the normalization was constrained to be <1 x 10"^^, resulting in 
a distance constrain of >316 Mpc. This unrealistic distance limit further demonstrates that 
a simple NSA model does not provide a good fit to the data. 

We also fitted the data using the NSA model but with the normalization fixed and 
assuming three distances: 5, 10, or 15 kpc. The results are listed in Table 1 and shown 
in Figure 3 (top). From this table it can be seen that such a model does not provide an 
acceptable fit to the data since 100% of the simulated spectra have lower values for the 
Cash-statistics than the data themselves. We also fitted a simple blackbody model to the 
data but again such a model did not provide an adequate fit to the data with a fit quality 
of 0.94 and a kT of O.Gtoi keV (and iVn < 3 x 10^° cm"^). The radius of the emitting area 
was < 0.12 X yq-y7 km, with d the distance in kpc, much smaller than the expected radius 
of a neutron star. The inability of the thermal models to provide an adequate fit to the 
data and the fact that four photons are detected with energies above 3 keV point to a X-ray 
spectrum which is significantly harder than a simple thermal model. Therefore, we fitted the 
X-ray spectrum with a power-law model (Fig. 3; Tab. 1). The fit quality of 0.68 using such a 
model is considerably better than for a thermal model (0.68 versus >0.88). The fit resulted 
in a power-law photon index F of 1.8!to'5, a column density of < 6 x 10^*^ cm~^, and a 0.5-10 
keV flux of 6^2 ^ 10~^^ ergs s^^ cm^^ (corrected for absorption). This gives an unabsorbed 
0.5-10 keV X-ray luminosity of 7^1 x 10^^( ^q'^ ^ )^ ergs s^^, with d the distance in kpc. 

To obtain limits on the temperature and luminosity of a possible thermal component in 
the data, we fltted the spectrum with a power-law component plus a neutron-star atmosphere 
component. We again flxed the normalization of the neutron-star atmosphere component 
to correspond to 5, 10, or 15 kpc but left the other parameters free (i.e., the column den- 
sity, the effective temperature, the photon index, and the normalization of the power-law 
component). The maximum allowed temperatures and bolometric luminosities are listed in 
Table 2, together with the minimum fraction of the 0.5-10 keV flux which is due to the 
power-law component. The table shows that the maximum allowed kT^o is between 0.04 and 
0.05 keV with a corresponding maximum allowed bolometric luminosity of 0.4 — 1.7 x 10^^ 
ergs s^^ (for both parameters the highest values are obtained for the cases which assume 
the largest distance). For all assumed distances, at least ~70% of the flux in the 0.5-10 keV 
range is due to the power-law component. This demonstrates that the power-law component 
dominates the X-ray emission in the range 0.5-10 keV for XTE J0929-314. 



2.1.4- Spectral variability 

In § 2.1.2 we demonstrated that XTE J0929-314 exhibited variabihty in its count rate 
at the 95% confidence leveL It is possible that this count rate variability is accompanied 
by or even due to changes in the X-ray spectrum of the source. If true, the spectral results 
we obtained in § 2.1.3 only represent a time-averaged view of the spectral properties of the 
source. To investigate if indeed XTE J0929-314 exhibited spectral variability we calculated 
X-ray hardness ratios from the first part of the observation (between and 10 ksec after the 
start of the observation; 13 photons were detected in this time interval) and the second part 
of the observation (the remaining ~14 ksec of the observation; 9 photons were detected). 
As X-ray hardness ratio we used the logarithm of the ratio between the number of photons 
detected with energies >1 keV and the number of photons detected with energies <1 keV. 
In Figure 4 [top panel) we show the hardness ratio versus time. During the second part of 
the observation the hardness ratio is smaller than during the first part, but the error bars 
are large and the data are consistent with a non-variable spectrum. 

Hong, Schlegel, & Grindlay (2004) have argued that the X-ray hardness ratios normally 
used in the literature (e.g., those we used in the previous paragraph) are not appropriate for 
sources with low numbers of counts. Since we have roughly 10 counts to construct our hard- 
ness ratios, these ratios might not represent an accurate description of the spectral shape of 
XTE J0929-314. Therefore, we also calculate so-called "quantiles" using the quantile method 
outlined by Hong et al. (2004)^. Hong et al. (2004) found that the quantity /ofi'io j^ A" , with 
Qso the median quantile, provided a good indication for the spectral hardness of a source. 
Therefore, in the bottom panel of Figure 4 we plot this quantity as a function of time. Sim- 
ilar to the hardness ratios, the error bars on the quantiles are such that a constant spectral 
shape throughout the observation cannot be excluded. 

As a last way of investigating possible spectral variability, we extracted the X-ray spec- 
trum for each of the two time intervals used above. The results of the spectral fitting are 
listed in Table 3. When fitting the two data sets separately, a power-law model produces 
acceptable fits (with fit quality of 0.6-0.7) with photon indices of 1.3^g;7 and 2.6lQ;g for the 
first and second part of the observation, respectively. The values again indicate that the 
source spectrum had softened during the course of the observation, although the indices 
are consistent with each other (within the errors). The 0.5-10 keV X-ray flux decreased 



^The errors on the quantiles were calculated using the method of Maritz & Jarret (1978; see Hong 
et al. 2004) but Hong et al. (2004) demonstrated, using simulations, that the errors on the quantiles are 
overestimated when the source counts are very low, as in our case. Therefore, we corrected the quantile 
errors using the corrections provided by Hong et al. (2004). 
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by a factor of 3-4 during the observation indeed suggesting variability in the brightness of 
the source (note, again the errors are large and only minimal variability might have been 
present). To further investigate how significant the softening of the X-ray spectrum is, we 
fitted the two data sets simultaneously. During this fit, we assumed that the column density 
did not vary during the observation and we tied it between the two data sets. The fit re- 
sults were consistent with those obtained when fitting the data sets separately (Tab. 3). In 
Figure 5 we plot the confidence regions of the photon indices, showing that the two indices 
are only different from each other at the <90% confidence level. This is illustrated in the 
figure by the line of equal index between two data sets which intersects the 90% confidence 
contour. We conclude that the limited statistics does not allow for any strong conclusions 
on the possible softening of the X-ray spectrum suggested by Figure 2. 

We investigated if a thermal model could fit the two data sets, but when using a NSA 
model with fixed normalization (to a distance of 5, 10, or 15 kpc), the fit qualities were 
always 1.00 for the first part of the observation and between 0.89 and 0.94 (larger for smaller 
assumed distance) for the second part of the observation. Also a simple black-body model 
does not produce a good fit with fit qualities of 0.85 and 0.87 for the first and second part of 
the observation, respectively. This clearly demonstrates that during the whole observation, 
the source had an X-ray spectral shape which was significantly different (i.e., harder) from a 
simple thermal model. This confirms our conclusions in § 2.1.3. Investigations into whether 
or not the fractional contribution of the power-law component to the 0.5-10 keV fiux changed 
during the observation did not result in useful insights due to the limitation on the statistics 
of the individual data sets. 



2.2. XTE J1751-305 

We begin a similar analysis for XTE J1751-305 by searching for point sources in our 
data (using 'wavdetect'). We detected 34 sources in our image but XTE J1751-305 itself was 
not detected (a discussion of the properties of these 34 sources, unrelated to XTE J175 1-305, 
is beyond the scope of this article). In Figure 6 we show the 0.5-7 keV X-ray image near 
the position of XTE J1751-305. Several point sources are detected close to the position of 
XTE Jl 75 1-305 (obtained during outburst using Chandra] Markwardt et al. 2002), but only 
1 photon was detected (in the energy range 0.5-7 keV) in a circle with a radius of 1" around 
the position of XTE J1751-305 (for Chandra^s point-spread-function, 1" corresponds to the 
encirclement of 90% of the energy), which does not constitute a significant detection. When 
using also the data at lower energies (i.e., using the energy range 0.3-7.0 keV), we still could 
only detect 2 photons within an 1" radius. When using a larger radius of 1.5", 3 photons 



-10- 



are detected for the energy range 0.5-7 keV, and 5 photons for the energy range 0.3-7 keV. 
However, in these cases more than half the possible source photons would fall outside the 
circle which should encircle >90% of the energy, which is unlikely. Moreover, the positions 
of the two softest photons (those with energies <0.45 keV) were consistent with the position 
of a bright K star close to the position of XTE J175 1-305 (1.2" away; the two photons have 
a distance of 0.5" and 0.8" from the position of this star). Thus these two photons could 
have possibly originated from this star and not from XTE J1751-305, especially because the 
high column density toward the pulsar (~10^^ cm~^; Miller et al. 2003) would make unlikely 
any detection of its photons which have energies below 0.5 keV. 

The closest, significantly detected X-ray source (the one to the south-east of XTE J1751- 
305) is located 3.8" away^ from XTE J1751-305. The absolute positional accuracy of Chandra 
is 0.6" (90% confidence level; 0.8" for 99% confidence level^) strongly suggesting that this 
south-east source is not the quiescent counterpart of XTE J1751-305. It might be possi- 
ble that during our observation the positional accuracy of Chandra was (for some reason) 
extremely inaccurate. To test this hypothesis, we compared our X-ray image with the K- 
band image of the region around XTE J1751-305 published by Jonker et al. (2003). Of 
the 13 X-ray sources located in the region covered by this X-band image, 8 have K objects 
located within a radius of 0.6" (including the X-ray source located 3.8" away from XTE 
J1751-305 as well as the four other reasonably bright X-ray point sources visible in Fig. 6). 
We simulated 10,000 X-ray point sources with random positions in the field of the K image. 
Only <7% of the simulated sources have a K star (of similar or higher brightness than the 
potential K counterparts of the detected X-ray point sources) falling within a 0.6" radius of 
their positions. The coincidence rate between the detected X-ray sources and the K stars is 
significantly higher (~60%) than the simulated one demonstrating that most of the K stars 
located in the X-ray error circles are indeed the infrared counterparts of those X-ray sources. 
Our analysis strongly suggests that the positional accuracy of our X-ray image is within the 
standard accuracy quoted for Chandra. Therefore, we conclude that XTE J1751-305 was 
not conclusively detected during our Chandra observation. 

To calculate a flux upper limit for XTE J1751-305, we assume that <5 photons^ (for 
the energy range 0.5-7 keV) were detected from the source, resulting in a count rate upper 



■^The coordinates of this source are: right ascension 17^ 51™ 13". 673, decHnation —30° 37' 26". 40. The 
coordinates are for epoch J2000.0 and have errors of 0.6" (90% confidence levels). 

^See http://cxc.harvard.edu/cal/ASPECT/celmon/ 

^Using 5 photons as an upper limit can be interpreted (Gchrels 1986) as a 97% or 90% confidence upper 
limit for 1 or 2 detected photons from the source, respectively, which is the number of photons we possible 
detected from XTE J 175 1-305. 
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limit of 1.16 xlO-'' counts s"^ (0.5-7 keV). We used PIMMS^ to calculate the flux upper 
limits assuming different spectral model. As column density we used the value as observed 
during outburst (A^h = 9.8 x 10^^ cm~^; Miller et al. 2003) and we assume a black-body 
spectral shape^ (with temperatures kT of 0.05-0.7 keV) or a power-law spectral shape (with 
a photon index F of 0-4). Using a black-body model, we find that the 0.5-10 keV flux upper 
limit Fupper (in units of ergs s^^ cm^^) can be accurately approximated by 

kT 

logioF^pp^r = -14.787 + 4.09e~o-o827 kcv (1) 

Using a power-law model, we flnd 

i^upper = 5.76 X 10"^^ - 3.05 X lO^^^r + 7.41 x lO^^^F^ (2) 

For a reasonable range in kT of 0.1-0.3 keV (as observed in other quiescent neutron star 
X-ray transients), we flnd 0.5-10 keV upper limits of 0.2 — 2.7 x 10"^^ ergs s^^ cm^^ (the 
highest upper limits are obtained when the assumed kT is smallest). For a reasonable F of 
1-3, the 0.5-10 keV flux upper limit is 2.6 — 3.5 x 10^^^ ergs s^^ cm^^ (the minimum upper 
limit is obtained when F ~2). 



3. Discussion 

We have observed the accretion-driven millisecond X-ray pulsars XTE J0929-314 and 
XTE J1751-305 in their quiescent states. This triples the number of such sources to be 
observed in their quiescent states (the flrst one being SAX J1808. 4-3658; Stella et al. 2000; 
Dotani, Asai, & Wijnands 2000; Campana et al. 2002). We could not detect XTE J1751-305 
during our observations, which resulted in upper limits on the quiescent X-ray flux (0.5-10 
keV; corrected for absorption) of this system of 2.6 — 3.5 x 10~^^ ergs s~^ cm~^ when we 
assume that the source exhibited a power-law spectral shape with photon index between 1 
and 3, or 0.2 — 2.7 x 10^^^ ergs s~^ cm~^ when assuming a black-body spectral shape with kT 
of 0.1-0.3 keV. We conclusively detected XTE J0929-314, but we detected only 22 photons 
(energy range 0.3-8 keV) during our observation which limits any detailed analysis. The 



® Available http://heasarc.gsfc.nasa.gov/Tools/w3piinms.htinl 

^Although a neutron-star atmosphere model is usually fitted to the X-ray spectra of quiescent neutron 
star X-ray transients, such a model cannot be used in PIMMS and we had to resort to the black-body model 
as an approximation. 
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X-ray spectrum could not be fitted with a simple thermal model alone (such as a black-body 
or a neutron-star atmosphere model) but a single power-law model provided an adequate 
fit to the data resulting in a 0.5-10 keV X-ray luminosity of 7^1 x 10^^( ^q'^ ^ )^ ergs s~^ 
(with d in kpc). The upper limits on the effective temperature of a thermal component 
ranged from 0.04 to 0.05 keV with corresponding upper limits on the bolometric luminosity 
of 0.4-1.7 xlO^^ ergs s"^ (both quantities are for an observer at infinity). We determined 
that at most ~30% of the 0.5-10 keV fiux could have come from such a thermal component. 
We could demonstrate that the source was variable at the 95% confidence level but, within 
the limitation of the statistics, we could not find evidence for spectral variability during 
our observation. Next, we will first discuss what can be inferred from the results of XTE 
J0929-314 before we discuss the non-detection of XTE J1751-305. 



3.1. The power-law component of XTE J0929-314 

Our results suggest that the X-ray emission of XTE J0929-314 in its quiescent state is 
dominated by a power-law component and not by a thermal component as is usually observed 
during the quiescent states of neutron-star X-ray transients. This makes XTE J0929-314 
similar to SAX J1808. 4-3658, which was until now the only accreting millisecond pulsar to 
have been observed in its quiescent state. In quiescence, SAX J1808. 4-3658 was found to 
have a 0.5-10 keV X-ray luminosity of ~5 x 10^^ ergs s^^ when fitted with a power-law 
model with index of 1.5 ± 0.3 (Campana et al. 2002). This is rather similar to the 0.5-10 
keV luminosity (~ 7 x 10^^ ergs s~^ for a distance of 10 kpc) and power-law index (1.8!']q'5) 
of XTE J0929-314. Furthermore, the quiescent X-ray emission between 0.5 and 10 keV 
was in both sources dominated by the power-law component (>90% in SAX J1808. 4-3658 
and >70% in XTE J0929-314). Although our statistics on the XTE J0929-314 data are 
poor, the only difference between the two sources seems to be that SAX J1808. 4-3658 has a 
somewhat lower luminosity in quiescence and perhaps a bit harder spectrum. We note that 
the errors on the parameters also suggest that both sources could have equal luminosities 
and spectral hardness. Furthermore, XTE J0929-314 exhibited luminosity variability during 
our observation possibly accompanied by (or even due to) spectral variability. This means 
that XTE J0929-314 could be at times intrinsically fainter than SAX J1808. 4-3658 or have a 
harder spectrum (note that it is less likely that XTE J0929-314 was fainter and harder than 
SAX J1808. 4-3658, since the indication is that the hardest spectrum was observed when 
XTE J0929-314 was brightest). 

The results on XTE J0929-314 and SAX J1808.4-3658 might suggest that a power- 
law dominated X-ray spectrum in quiescence and a low quiescent 0.5-10 keV luminosity 
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are common properties of accretion-driven millisecond X-ray pulsars. However, recent work 
(as discussed below) on two other weak quiescent neutron-star X-ray transients, which do 
not exhibit pulsations, suggests that this might instead be a property of many weak qui- 
escent neutron-star X-ray transients in general and not only of the accreting millisecond 
X-ray pulsars. Jonker et al. (2004a) and Tomsick et al. (2004) reported on quiescent X-ray 
observations of SAX J1810. 8-2609 and XTE J2123-058, respectively. For both sources it 
was found that their quiescent X-ray spectra could be fitted with a power-law model with 
indices of ~3 and an 0.5-10 keV X-ray luminosity around ~10^^ ergs s~^. No thermal com- 
ponent could be detected with a maximum contribution to the 0.5-10 keV flux of ~50% 
in SAX J1810.8-2609 (Jonker et al. 2004a) and ~60%-70% in XTE J2123-058 (Tomsick 
et al. 2004). These luminosities are somewhat higher than those of SAX J1808. 4-3658 and 
XTE J0929-314 and the X-ray spectra somewhat softer, but they are among the lowest 
luminosity quiescent neutron-star X-ray transients for which no X-ray pulsations were seen 
during their outbursts. 

Jonker et al. (2004a, b) used the results obtained for SAX J1808.4-3658, SAX J1810.8- 
2609, and XTE J2123-058 together with those obtained for the other quiescent neutron star 
X-ray transients, to show evidence that, in their quiescent spectra, the contributions of the 
power-law components to the 0.5-10 keV flux seem to increase (from at most several tens 
of percents to over 90%) when the quiescent source luminosities decrease from ~10^^ ergs 
s^^ to approximately 5 x 10^^ ergs s^^. They also found a similar trend for higher quiescent 
luminosities of the transients: the fractional power-law contribution also increases when the 
quiescent source luminosities increase from the comparison luminosity of 1 — 2 x 10^^ ergs s^^. 
Our results show that XTE J0929-314 fits between SAX J1808.4-3658 on the low-luminosity 
side, and SAX J1810. 8-2609 and XTE J2123-058 on the higher-luminosity side, which pro- 
vides additional support to the possible correlations found by Jonker et al. (2004a, b). This 
result would suggest that hard quiescent emission is not a feature unique to accretion-driven 
millisecond X-ray pulsars but rather a general feature of quiescent neutron-star X-ray tran- 
sients when they are at low quiescent luminosities (with increasingly harder spectra when 
the luminosities get lower). The current data still suggest that the accretion-driven millisec- 
ond X-ray pulsars might have slightly lower luminosities and harder spectra than the other 
sources, but it remains to be seen if this holds when additional neutron-star X-ray transients 
are observed in their quiescent states. Additional neutron star X-ray transients need to be 
studied in their quiescent states to better understand the power-law spectral component in 
the quiescent spectra. Such studies might also clarify if the quiescent transient in Terzan 
5 (EXO 1745-248; Wijnands et al. 2004), which does not follow this correlation (Jonker et 
al. 2004b), is unique among the quiescent neutron star systems, or if it is the first example 
of many similar systems which would suggest that this correlation is spurious. 
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3.2. The non-detection of a thermal component in XTE J0929— 314 

We can investigate whether our upper limits on the contribution of the thermal com- 
ponent to the quiescent spectrum of XTE J0929-314 are consistent with what would be 
expected from the cooling neutron star model proposed by Brown et al. (1998). The quies- 
cent thermal luminosity predicted by this model depends on the time-averaged (over >10,000 
years) accretion rate of the source. For XTE J0929-314 we can estimate its time-averaged 
accretion rate in two ways; either by using the outburst fluence of the source and the limited 
information about the outburst recurrence time, or by assuming that the mass transfer is 
solely due to gravitational radiation and all matter is eventually accreted onto the neutron 
star. 



3.2.1. Constraints obtained when using the 2002 outburst fluence 

According to Galloway et al. (2002), the 2-60 keV fluence of the 2002 outburst of XTE 
J0929-314 was 4.2 x 10"'^ ergs cm^^, with an estimated bolometric correction of ~2.34. This 
results in a bolometric fluence of ~9.8 x 10~^ ergs cm~^. During the lifetime of RXTE 
(which is at present 8.5 years), only 1 outburst has been observed, resulting in a time 
averaged accretion flux (Face) of < 3.7 x 10~^^ ergs s~^ cm~^. Using the Brown et al. (1998) 
model and assuming standard core cooling, the expected quiescent flux Fq is approximately 
given by Fq = ^^ (Brown et al. 1998; Wijnands et al. 2001; Rutledge et al. 2002). For 
XTE J0929-314, this results in Fq < 2.7 x 10^^^ ergs s^^ cm^^, which is consistent with the 
upper limits on the bolometric fluxes of a possible neutron-star atmosphere component in 
the quiescent spectra (which are <1.5 xlO^^^ ergs s^^ cm^^). However, in this scenario the 
source would have to remain quiescent for a signiflcant period (i.e., much longer than the 
8.5 years we assumed here) to allow the predicted limits to come down to the measured flux 
limits. 

Using our measured flux limits, we can actually estimate how long XTE J0929-314 
has to be in quiescence for the neutron star surface to be as cool as we measure if the 
neutron star cools down using standard core cooling processes. From Fq = -^-^^ it can be 
derived (Wijnands et al. 2001) that Fq ^ ^-^ x ■^, with (Fq) the averaged flux during 
outburst, to the averaged time the source is in outburst, and tq the averaged time the source 
is in quiescence. From Galloway et al. (2002), it follows that to ~ 73 days (if the 2002 
outburst had a duration which is typical for the source) and (Fq) = 1.6 x 10^^ ergs s^^ 
cm~^. Using our limits on Fq < 1.5 x 10^^"^ ergs s~^ cm~^, then tq > 5.8 x 10^ days, 
resulting in tq > 160 years. Although this predicted quiescent period is quite long compared 
to the known recurrence time of several other neutron-star X-ray transients (ranging from 
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less than 1 year up to several decades; e.g., Chen, Shrader, & Livio 1997), the disk instabihty 
model proposed for the outburst behavior of X-ray transients (Lasota 2001) might be able 
to accommodate such long quiescent episodes. We note that due to the very low peak fluxes 
during the 2002 outburst, we cannot rule out that, prior to the launch of RXTE, one or more 
outbursts occurred during the last four decades (i.e., since the birth of X-ray astronomy). 
Weak outbursts, such as the 2002 outburst of XTE J0929-314, could have easily been missed 
by X-ray satellites other than RXTE. 



3.2.2. Constraints obtained when assuming mass transfer solely via gravitational radiation 

Galloway et al. (2002) stated that if the mass transfer in XTE J0929-314 is driven by 
gravitational radiation then the time-averaged mass transfer rate Mgr in this system should 
be 

Mgr = 5.5 X 10-12(^^)1(^^)2 Mo yr'^ (3) 

h.4 Mq^ ^0.01 Mq^ ^ ^ ^ ^ 

with Mns the mass of the neutron star and M^ the mass of the companion star (both in solar 
masses). From the model by Brown et al. (1998; see also Rutledge et al. 2002) the quiescent 
luminosity Lq can be estimated via 

with (M) the time-averaged mass accretion rate onto the neutron star, and Q the amount 
of heat deposited in the crust per accreted nucleon (e.g., Haensel & Zdunik 1990). If we 
assume that the matter transfer from the companion star is solely driven by gravitational 
radiation and that all the matter transferred from the companion star is eventually accreted 
onto the neutron star then (M) = Mqr. Substituting equation 3 in equation 4 gives 

Lq = 4.8 X 1032( J^)l(^^)2 ^ ergs s'^ (5) 

^ ^I.AMq' ^O.OIMq' 1.45 Me V ^ ^ ' 

If Mjis = 1.4 Mq, then Galloway et al. (2002) obtained a minimum mass for the companion 
star of 0.008 Mq (thus M^ > 0.008 Mq). If we further assume Q = 1.45 MeV, then 
Lq > 3 X 10^2 ergs s^^. Only for large distances (>15 kpc) can the upper-limits on any 
contribution of a thermal component be consistent with this lower-limit on the quiescent 
luminosity (if M^ is not much greater than the minimum allowed value of 0.008 M0). For 
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smaller distances the predicted quiescent luminosity for a thermal component is a factor 2-8 
times too high (this discrepancy becomes worse if M^ becomes larger). Firm conclusions can 
only be obtained when the exact luminosity of the thermal component can be measured with 
better data. However, our results already suggest that the quiescent luminosity predicted is 
higher than that measured if we assume that the matter transfer from the companion star is 
solely driven by gravitational radiation and that all the matter transferred is accreted onto 
the neutron star. Possible reasons for this discrepancy are that not all the matter transferred 
is accreted onto the neutron star or that enhanced cooling processes occur in the neutron-star 
core. 



3.3. The non-detection of XTE J1751-305 

Despite the fact that we could not conclusively detect XTE J1751-305 during our obser- 
vation, the limits on the X-ray flux of this source can still be used to provide further insight 
into the quiescent emission of neutron star X-ray transients. We do not know the spectral 
shape of the quiescent emission of this system and therefore in § 2.2 we calculated the flux 
upper limits assuming two different spectral models (either a black-body or a power-law 
model; equations 1 and 2). The limits obtained using the black-body model are most useful 
for constraining the cooling neutron star model which assumes that the emission is due to 
thermal radiation from the neutron-star surface. We will attempt to constrain this model 
once again taking into account the limited observational knowledge about the time-averaged 
accretion rate of the source (§ 3.3.1) and assuming that the mass transfer is solely due to 
gravitational radiation and that all matter will eventually be accreted by the neutron star 
(§ 3.3.2) . Here we first discuss the possibility that the source emitted a power-law shaped 
quiescent X-ray spectrum. 

The two other accretion-driven millisecond X-ray pulsar which have been observed in 
quiescence have shown a power-law dominated quiescent X-ray spectrum with photon indices 
of 1.5±0.3 (SAX J1808.4-3658; Campana et al. 2002) and 1.8±0.6 (XTE J0929-314; Tab. 1). 
If we assume that XTE J1751-305 has also a power-law shaped quiescent spectrum, then, 
using equation 2, we get 0.5-10 keV flux upper limits on the quiescent emission of 2.5 — 
3.5 X 10^^^ ergs s^^ cm^^ (assuming a range of 1-3 for the photon index, which covers the 
allowed values of the photon indices for SAX J1808.4-3658 and XTE J0929-314), giving 0.5- 
10 keV luminosity upper limits of 1.9 — 2.7 x 10^^( g ^ ^ )^ ergs s~^ (d in kpc). If the distance 
toward XTE J1751-305 is <11 kpc (thus the source is located in the Galactic center region) 
and the quiescent spectrum follows a power-law shape, then XTE J 175 1-305 would be the 
intrinsically faintest quiescent neutron-star X-ray transients so far known. However, if the 
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source is located at a larger distance, SAX J1808. 4-3658 would still remain the faintest 
quiescent system (with a luminosity of ~ 5 x 10^^ ergs s^^; Campana et al. 2002), unless it 
can be confirmed that XTE J0929-314 can become at times fainter than SAX J1808. 4-3658 
(see 5 3.1). 



3.3.1. Constraints obtained when using the 2002 outburst fluence 

According to Markwardt et al. (2002), the 2-200 keV fluence of the 2002 outburst 
of XTE J1751-305 was 2.5 x 10~^ ergs cm~^. Since no bolometric correction factor is 
currently available, we will use this fluence, but we note that the bolometric fluence could 
be significantly higher (a factor of a few) and thus also the predicted quiescent thermal flux. 
Markwardt et al. (2002) found a previous outburst of the source in June 1998, resulting in 
a recurrence time of ~3.8 years. If this is the typical recurrence time of the source, then the 
time-averaged accretion rate will be (-Face) ~ 2.1 x 10^^^ ergs s^^ cm^^. Since the predicted 
quiescent flux is given by Fq = -^^ (§ 3.2.1), this results in Fq = 1.6 x 10^-*^'^ erg s~^ 
cm~^, which is significantly larger than the maximum bolometric flux (2.7 x 10"^"^ ergs s~^ 
cm~^) we obtained for a possible thermally shaped quiescent spectrum. It might be possible 
that the average recurrence time of the source is larger than the 3.8 years found for the two 
outbursts so far observed. Using the method outlined in § 3.2.1, we find that the source 
must have typical quiescent episodes of >22 years for the predicted quiescent flux to come 
down to the calculated flux limits. Such a recurrence time scale is not uncommon among the 
neutron-star X-ray transients. Another possibility for the discrepancy between the predicted 
quiescent flux and the limits on the actual flux is that the neutron star in XTE J1751-305 
is colder than expected due to enhanced cooling in its core. 



3.3.2. Constrains obtained when assuming mass transfer solely via gravitational radiation 

Similar to XTE J0929-314, it also has been suggested for XTE J1751-305 (Markwardt 
et al. 2002) that the mass transfer might be due to gravitational radiation. If true, then 
Mqr in this system should be (Markwardt et al. 2002) 

Mgr = 1.2 X 10-"(^^^)i( ^ f Mr, YT-\ (6) 

^^ ^1.4 Mq' ^0.0137 Mq' ^' 

If we assume that the matter transfer is solely driven by gravitational radiation and that all 
the matter is eventually accreted onto the neutron star then we can again set (M) = Mgr. 
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Substituting equation 6 in equation 4 gives for the predicted quiescent luminosity (again 
assuming Q = 1.45 MeV) 

Lq = 1.0 X 1033( J^)l( M^ )2 ergs s'^. (7) 

Markwardt et al (2002) obtained a minimum mass for the companion star of M^ > 0.0137 Mq, 
when assuming a neutron star mass of 1.4 Mq. Using these values for the masses of the com- 
panion and the neutron star, we find that the predicted quiescent luminosity is given by 
Lq > 1 X 10^^ ergs s^^. The 0.5-10 keV flux upper limit assuming a black body spectral 
model is given by equation 1. If kT > 0.1 keV (which is typically observed for other neutron 
star X-ray transients in their quiescent states), this gives a 0.5-10 keV luminosity upper 
limit of > 2 X 10^^ s "^k c ^^Ss s~^ {d in kpc). The bolometric luminosity limit could be a 
factor of a few larger. Only for relatively large distances (> 10 — 15 kpc; depending on the 
bolometric correction factor) can the measured upper limit on the luminosity become consis- 
tent with the predicted luminosity. For effective temperatures below 0.1 keV, the bolometric 
luminosity as calculated via Lboi = 47rcr i?^Tjj (with R and T the neutron star radius and 
surface temperature for an observer at infinity and a the Stefan-Boltzmann constant) will 
become quite low and will not be consistent with the predicted luminosity. This indicates 
that the effective temperature cannot be too low in order to be consistent with the predicted 
quiescent luminosity. 



3.4. Constraining the magnetic field of the neutron stars in XTE J0929— 314 

and XTE J 175 1-305 

Recently, Burderi et al. (2002) and Di Salvo & Burderi (2003) constrained the magnetic 
field strengths of the neutron stars in several neutron-star X-ray transients (i.e., KS 1731- 
260, SAX J1808. 4-3658 and Aql X-1) based on their measured quiescent luminosities and 
on the knowledge of the spin rate of their neutron stars^. To constrain the magnetic field 
strength, they used two proposed mechanisms which might produce X-rays in quiescence: 
a) residual accretion onto the neutron star at a very low accretion rate; b) neutron star 
rotational energy which is converted into radiation (a fraction of which might be released 



^SAX J1808. 4-3658 exhibits coherent osciUations in its persistent X-ray emission during outbursts which 
reflect the spin rate of the neutron star (Wijnands & van der KUs 1998). KS 1731-260 and Aql X-1 exhibit 
nearly-coherent oscillations during thermonuclear flashes on the neutron star surfaces (Smith, Morgan, & 
Bradt 1997; Zhang et al. 1998). Those oscillations are likely to be directly related to the neutron star spin 
rate as well (see, e.g., Strohmayer & Bildsten 2004 for a recent review). 
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in X-rays) because of the presence of a rotating magnetic dipole. They note that thermal 
emission from a coohng neutron star which is heated during outburst might also contribute 
to the quiescent X-ray emission and that the measured quiescent luminosities are therefore 
an upper limit on the X-ray contribution due to the above discussed scenarios. 

They divided scenario a into two possibilities. The first (al) assumes that the radius 
of the magnetosphere is inside the co-rotation radius and accretion onto the neutron star 
surface is possible; the second {a2) assumes that the magneto-spheric radius is outside the 
co-rotation radius but inside the radius of the light cylinder and accretion onto the neutron 
star is not possible but an X-ray emitting accretion disk might still exist with an inner radius 
larger than the co-rotation radius. Scenario b can also be divided into two possibilities. The 
first (bl) assumes that the X-ray emission in quiescence is due to reprocessing of part of 
the bolometric luminosity from the rotating neutron star into X-rays in the shock front 
between the pulsar wind and the circumstellar matter; the second {b2) assumes that the 
X-ray emission is intrinsic emission from the radio pulsar. 

Following the method outlined by Burderi et al. (2002) and Di Salvo & Burderi (2003) 
and using the maximum quiescent luminosity from our fits (1.2 x 10^'^{ "^^ ^ )^ ergs s~^, with 
d in kpc) and the spin rate (185 Hz; Galloway et al. 2002) of XTE J0929-314, we found that 
the magnetic field strength of the neutron star in this system should be < 2 x 10'^ ^^'^ Gauss 

(scenario al), < 3x lO^yofc ^^^^^ («^)' < ^^ lO^iofc ^^^^^ (^-^)' ^^^ < 3x 109(^q^)0-™ 
Gauss {b2). Therefore, from these four scenarios the magnetic field strength of the neutron 
star in XTE J0929-314 can be constrained to be < 3 x IO^tptt — Gauss. 

10 kpc 

We also calculated the limits on the magnetic field strength of the neutron star in XTE 
J1751-305. The maximum fiux limit of 2.7 x 10~^^ ergs s~^ cm~^ we obtained was for a 
spectral model which we assumed to be a simple black body shape, resulting in a luminosity 
limit of 2 X 10^^ g ^ ^ ergs s~^ (with d in kpc). Using the spin frequency of 435 Hz (Markwardt 
et al. 2002), this gives limits on the magnetic field strength of < 1 x lO'^ g |^ ^ Gauss (scenario 
al), < 5 X 10^^ Gauss (a2), < 1 x 10^^ Gauss (61), and 7 x 10^^ Gauss (62). 
When assuming a power-law spectral model for the quiescent spectrum of XTE J1751-305, 
the upper limits on the fiuxes are significantly smaller, with a maximum luminosity limit 
of 3 X 10^^ g ^ c (with d in kpc), which would make the magnetic field limits a factor of >2 
less (the maximum upper limit would be again for scenario b2 and is < 3 x lO'^ g j^ ^ Gauss). 
Therefore, from these four scenarios the magnetic field strength of the neutron star in XTE 
J1751-305 can be constrained to be < 3 — 7 x lO^g-^-^ Gauss (depending on assumed spectral 
model for the quiescent spectrum). 

Although this method still has significant uncertainties (i.e., due to the lack of under- 
standing of the structure of the accretion fiow geometry in quiescent X-ray transients; see 
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Burderi et al. 2002 and Di Salvo & Burderi 2003 for a more detailed discussion), our limits 
on the magnetic field of the neutron stars in XTE J0929-314 and XTE J1751-305 are con- 
sistent with the expectation that these neutron stars have a low but non-negligible (since 
they exhibit pulsations during accretion outbursts) magnetic field strength. 

We thank Peter Jonker for allowing us to use his K-hand image of the region around 
XTE J1751-305 and for comments on an early version of this paper. We thank Jae Sub 
Hong for useful discussion about the quantile analysis method. 
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Fig. 1. — The Chandra image of the field around XTE J0929-314 (left) and a close-up of 
XTE J0929-314 (right) with the radio error circle for the source obtained during outburst 
(Rupen et al. 2002). In both figures, north is up and east is to the left. In the left panel the 
arrow indicates which source is XTE J0929-314. The source located to the north-west of 
XTE J0929-314 is CXOU J092919. 1-312244 and the source to the south-east (in the bottom 
left corner) is CXOU J092924.4-312420. 
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Fig. 2. — The energy of the detected photons versus the time of detection (since the start of 
the observation). 
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Fig. 3. — The time-averaged spectrum of XTE J0929-314. In the top panel, we show the 
best power-law fit to the data (solid line) and the best NSA fit through the data with a fixed 
normalization (to a distance of 10 kpc; dashed line). In the bottom panel, we again show 
the best power-law fit (solid line) but now the best NSA fit with the normalization left free 
(dashed line). For display purposes, the data were rebinned so that each bin has 3 counts 
but the spectra were fitted without any rebinning. 
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Fig. 4. — The hardness ratio (top) and the log 



10 



i-c 



{bottom) versus the time during the 



observation. The hardness ratio is defined as the logarithm of the ratio between the number 
of photons detected with energies >1 keV and the number detected with energies <1 keV, 
and Qso is the median quantile as defined by Hong et al. (2004). 
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Fig. 5. — A comparison of the photon index of the first part of the data (0-10 ksec; Fi) with 
that of the second part of the data (the remaining ~14 ksec of the observation; F2). The 
contours are the 68%, 90% and 99% confidence levels. The cross marks the position of the 
best-fit values of the indices. The dashed line indicates where both indices are equal. 
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Fig. 6. — The 0.5-7 keV image of the region around XTE J1751-305. East is to the left 
and north is up. Several point sources are detected close to the position of XTE J1751-305 
but none of these sources is consistent with the position of XTE J1751-305 (as indicated by 
the circle which has a radius of 1", which means that 90% of the energy of the photons is 
encircled in this region). 
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Table 1. Spectral fits of XTE J0929-314 



Model" 



iVn 


kT^ 


r 






pb 






Fit 


quality' 


1022 cm-2) 


(keV) 




(10- 


-15 


ergs s" 


'^ cm" 


-') 






0.3i°:? 


0.045±0.005 


— 






9tl 








1.00 


0.2 ±0.2 


0.058t',Z 


— 






9tl 








1.00 


0.2l°i 


0.067±0.08 


— 






9^-i 








1.00 


<0.06 


— 


1 8+0-6 






6+^ 








0.68 



NSA {d = 5 kpc) 
NSA {d = 10 kpc) 
NSA {d = 15 kpc) 
Power-law 



" For the neutron-star hydrogen atmosphere (NSA) model we used that of Zavlin et al. (1996) 
with a mass of 1.4 Mq, a radius of 10 km, and a zero magnetic field strength for the neutron star 
in XTE J0929-314. The errors on the fit parameters are for 90% confidence levels. 

'' The 0.5-10 keV unabsorbed fiux. 

'^ The fit quality represents the fraction of 10,000 simulated spectra which have Cash-statistics 
better than that obtained for the actual data. 



-29- 



5 


<0.14 


1 O+0.9 


<0.04 


<0.4 


10 


<0.10 


1 1+^-° 
J--J--1.2 


<0.05 


<1.2 


15 


<0.08 


1 1+^-° 
J--J--1.3 


<0.05 


<1.7 



Table 2. Limits on any NSA contribution in XTE J0929-314" 

Distance d N-^ T kT^ Bolometric luminosity^ Fraction in the 

(kpc) (10^^ cm~^) (keV) (10^^ ergs s~^) power-law componenf^ 

>0.77 
>0.70 
>0.78 

" Results of spectral fits assuming a NSA plus a power-law model. For the NSA model we 
used that of Zavlin et al. (1996) with a mass of 1.4 M©, a radius of 10 km, and a zero magnetic 
field strength for the neutron star in XTE J0929-314. The upper limits are for 90% confidence 
levels. 

^ Bolometric luminosity of the neutron-star atmosphere component as seen by an observer at 
infinity 

■^ The minimum fraction of the 0.5-10 keV fiux in the power-law component 
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Table 3. Possible spectral variability in XTE J0929-314 



Time seled 
(ksec since 


;ion 

start) 


TVh 
(10^2 cm-2) 


r 

(10- 


-15 


Flux'' 

ergs s^^ cm"2) 


Fit quality 


Separate fits 


0-10 
10 - end 




<0.2 
<0.1 


1 O+0.9 

2 fi+^-2 




14+1^ 

3tl 


0.70 
0.62 


Simultaneous fits with tied A^h 


0-10 
10 - end 




<0.07 


1.3 ±0.7 

^•'-'-0.9 




14^- 

3tl 


0.66 



"■ The 0.5-10 keV unabsorbed flux 



